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Abstract 

Broad-band spectra of accreting X-ray binary pulsars can be fit by a phenomeno- 
logical model composed by a power law with a high energy rollover above 10 keV, 
plus a blackbody component with a temperature of few hundred eV. While, at least 
qualitatively, the hard tail can be explained in terms of (inverse) Compton scatter- 
ing, the origin of the soft component cannot find a unique explanation. Recently a 
qualitative picture able to explain the overall broad band-spectrum of luminous X- 
ray pulsars was carried out by taking into account the effect of bulk Comptonization 
in the accretion column. After a review on these recent theoretical developments, 
I will present a case study of how different modelization of the continuum affect 
broad features, in particular the cyclotron resonance features in Vela X-l. 
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1 Introduction 



Since the observation of pulsed emission fro m Centaurus X-3 (IGiacconi et al. 
1971 ) and the discovery of its binary nature ( Schreirer et al. . 1972f ). more than 



100 other X-ray binary pulsars (XBPs) have been observed. The qualitative 
physical s cenario able to ex plain their pulsed X-ray emission was first elab- 
orated by IShklovskiil ([19671 ) before their discovery: X-rays are produced in 
the conversion into radiation of the kinetic energy of matter coming from the 
stellar companion and accreting onto the neutron star (NS). Because of the 
interaction with the NS strong magnetic field, of the order of 10 n -10 13 G, 
matter is driven onto the magnetic polar caps, where it is decelerated and 
X-rays are produced. If the magnetic field axis is not aligned with the spin 
axis, the NS acts as a "lighthouse", giving rise to pulsed emission when the 
beam crosses our line of sight. For a detailed description of the XBP spectral 
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Fig. 1. Schematic block diagram of the physical processes occurring in a X-ray 
binary pulsar accreting fro m its companion . For a detailed discussion on the physics 
involved in each block see lOrlandinil (|2004h . 



properties it is therefore necessary to describe the interactions of the X-rays 
produced at the NS sur face with the hi ghly magnetized plasma forming the 
magnetosphere (see, e.g. Meszarosl . 19921 ) . This is a formidable task because we 
cannot use a linearized theory for the radiative transfer equations but we have 
to deal with the fully magnetohydrodynamical system. In Fig. 1 we present 
a very schematic block diagram of the physical scen ario that describ es the 
production and emission of pulsed X-rays from XBPs ([Orlandinil . 120041 ) . 



From the picture given above it should be clear that in order to be able to 
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describe the XBP spectra it is first necessary to understand how radiation 
is formed and then to study how this radiation interacts with the strongly 
magnetized plasma of the NS. As a consequence of this interaction, cyclotron 
res onance features (CRFs) are observed in XRB spectra (for a recent review 



sec 



Harding, 2003[ ). Because of their broad character, CRFs are quite sensible 



to the continuum modelization. After a review on the mechanisms of spectral 
formation in XBPs, I will present the case study of the Vela X-l CRFs and 
their dependence on the continuum adopted. 



2 Spectral formation mechanism in XBPs 



Once matter has penetrated the NS magnetosphere, it will follow the magnetic 
field lines up to the magnetic polar cap, where it will be decelerated (see 
Fig. 1). If the amount of matter falling onto the polar caps is high enough 
that an X-ray luminosity greater than about 10 37 erg/s is reached, then a 
radiation-dominated (collisionless) shock wil l form at some distance a bove 
the NS surface, creating an accretion column ( Basko and Sunvaev . 1976j ). For 



luminosities lower than 10 37 erg/s the radiation pressure is not sufficient to 
stop the infalling matter, which therefore can impact directly onto the NS. In 
this case the source of X-ray radiation is an accretion slab and emission occurs 
in a direction parallel to the magnetic field (pencil beam), while in the case of 
emission from an accretion column X-rays are emitted from the sides of the 
column (fan beam) because the column is optically thick to X-ray radiation. 
It is clear that the formation or not of an accretion column (that is, of a 
collisionless shock) depends on how much X-ray radiation is produced which, 
in turn, depends on the radiative transport through the infalling matter. In 
other words, the flow dynamics and the radiative transport are coupled, in a 
sort of "feedback effect". 

From a theoretical point of view, the main physical mechanism of emission 
in XBPs is (inverse) Compton scattering. According to the the value of the 
Comptonization parameter y, we expect to observe a (modified) blackbody 
spectrum if y <C 1, while for y ^> 1 inverse Compton scattering can be im- 
portant. If we define a frequency u co such that y{oj co ) = 1, then for uj uj cq 
the inverse Compton scattering is saturated and the emergent spectrum will 
show a Wien hump, due to low-e nergy photons up-scattered up to fvuj ~ 
3kT ( Rvbicki and Lightman . 1975f ). In the case in which there is not satu- 



ration a detailed analysis of the Kompaneets equation shows that the spec- 
trum will have the form of a power law modified by a hi gh energy cutoff 
( Rvbicki and Lightman . 19751 : Sunvaev and Titarchuk . 1980l ). 



Numerous attempts were m ade to numerically simulate XBP spectra (see e.g. 
Meszaros and Nagel 1985al fbh. While there was a qualitatively agreement with 
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the observations, an ad hoc source of soft photons was required in order to 
fit the low energy part of t he spectrum. This low energy thermal emission 
was discussed into detail bv lHickox et al. ( 2004 ). They find that, in the case 



of luminous (> 10 38 erg/s) pulsars, the most likely origin of this soft excess 
is the reprocessing of hard X-rays (produced by the NS) by optically thick 
material. This process can be excluded for less luminous sources (< 10 36 erg/s), 
for which the most likely origin of the soft emission is from the NS surface or 
circumstellar matter heated by the NS emission. For intermediate luminosities 
XBPs are likely present both kind of emission. 



Recently, iBecker and Wolfa (|2005aU bl) renewed the interest on spectral forma- 
tion in XBPs by introducing, in the presence of an accretion column, the effect 
of bulk or dynamical Comptonization. In the accretion column, the infalling 
electrons that act as scattering medium for the soft photons possess a pre- 
ferred motion: this implies that photons gain energy through first-order Fermi 
acceleration. In the case considered in the past of the thermal Comptonization 
process, photons gain energy via second-order Fermi acceleration because of 
the incoherent, stochastic motion of the plasma. Therefore the energization 
process occurs at the expense of the bulk kinetic energy of the infalling mat- 
ter and is not supplied by the gas internal energy. In this scenario the soft 
component comes out naturally from soft photons produced by the "thermal 
mound" at the base of the accretion column that are able to escape without 
experiencing many scatterings. 

Taking into account only bulk Comptonization and neglecting cyclotron emis- 
sion and absorption, Becker and Wolfi ( 2005bf ) were able to reproduce quali- 



tatively the general form of a XBP broad-band spectrum: they predict a high 
energy power law and a turnover at low energy. This model is not able to 
explain the high energy cutoff observed in many XBPs, but this is proba- 
bly due to the not inclusion of cyclotron processes: indeed the cutoff energy 
and the cyclotron resonance f eature (CRF) energy are strongly correlated 
( Makishima and Miharai Il992h . 



3 Comparison with observations 



From an observational point of view the phenomenological model able to de- 
scribe the broad-band spectra of accreting XBPs contains (i) a black-body 
component with temperature of few hundreds eV; (ii) a power law of photon 
index ~1 up to ~10 keV; and a (iii) a high energy (> 1 keV) cutoff that makes 
the spectrum rapidly drop above ~40-50 keV (see, e .g. lOrlandini and Dal FiumeL 
20011 : ICoburn et ail . 120021 iFilippova etHi l2005h . Superimposed there are 
emission line features due to fluorescence from ions at different ionization 
levels, and broad absorption "lines" due to cyclotron resonances. In Fig. 2 
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BeppoSAX SVP Observation of 4U 1626-67 
model: constant Tabs gausabs highecut (bbodyrad + powerlaw + gaussian) 
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Fig. 2. 0.1-100 keV pulse averaged spe ctrum of the XBP 4U 1626-67 as observed 



by BeppoSAX (|Orlandini et al.fl l998b). All the typical XBP spectral components 



described in the text are present in this source. 

we show the BeppoSAX 0.1-100 ke V, pulse averaged spectrum of the XBP 



4U 1626-67 (jOrlandini et a1.Ul998bh . in which all the components mentioned 



above are clearly present. 



The first descriptio n of the high energ y spectral rollover was quite crude 
( White et ah . 19831 ). and Tanaka (1986) "smoothed" the cutoff by introduc- 
ing the so-called Fe rmi-Dirac cuto ff. By analyzing a sample of XBP spectra 
observed by Ginaa. iMihara ( 199fih introduced the so-called NPEX (Negative 
Positive Exponential) model: 



NPEX(E) = (AE~ a + BE +P ) exp 



E 
~kT 



the components of which have a physical meaning becaus e, if = 2, it mimics 

the sa turated inverse Compton spectrum (see discussion in Orlandini and Dal Fiumel 
20011 ). Furthermore, because the (non relativistic) ener gy variation of a photon 
during Compton scattering is AE/E = {AkT—E)/mc 2 ( Rybicki and Lightmanl . 
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1975j ) then when E = E cyc , with E cyc the CRF energy, the medium is opti- 
cally thick and therefore E cyc ~ 4kT. In other words, the measurement of the 
CRF energy gives an order of magnitude estimate of the temperature of the 
electrons responsible of the resonance. 



3.1 The CRFs m Vela X-l 



The knowledge of the form of the continuum is of paramount importance for 
the determination of the physical parameters derived from broad features, as 
CRFs. In Fig. 3 we show the fits with different spectral models to the Bep- 
poSAX 0.1-100 keV, pulse-averaged spectrum of one of the best studied XBPs, 
Vela X-l. This source presents a controversial CRF at ~ 25 keV and a w ell 



deter mined CRF at ~ 50 keV ([Orlandini et all Il998at iKrevkenbohm et al. 
In order to check if the first CRF is real or an artifact due to a incor- 



rect modelization of the continuum, we fit the spectrum with three different 
continuum models: a broken power law (top panel in Fig. 3), and two cutoff 
power laws (second and third panel). In the first case we needed two CRFs, 
while in the other two cases only one CRF at ~ 50 keV was require d. In the 
secon d and third panel we modeled t he CRF with the CY CABS model (jMihara . 
1991), and a Gaussian in absorption (|Soong et allll99of ). respectively. In order 



to extract the CRF profile from the residual panel, we put to zero the CRF 
normalization. It is quite evident from the residual panels in Fig. 3 that the 
CRFs strongly depend on the choice of the continuum model (from an F-test 
they result statistically equivalent): the ~ 25 keV features is surely present 
with a broken power law continuum (although its energy is quite close to the 
break energy) while it disappears if using two cutoff power laws. In this latter 
case there is a hint of its presence by using a cycabs modelization of the 
CRF, while it is not present at all by using a Gaussian in absorption. A pulse 
phase spect roscopy performed on a longer BeppoSAX observation confirmed 
this result ( La Barbera et al. . 2003f ). We think that this feature could be due 
to a two-step steepening of the spectrum modeled with an i ncorrect smooth 
rollov er, as already observed in the XBP OAO 1657-415 f Or landini et al 
19991) . 



4 Conclusion 



After more that 30 years from their discovery, there is not yet a satisfactory 
theoretical model for the spectral emission in XBPs able to explain the ob- 
served spectra. Even if the recognition of the importance of bulk Comptoniza- 
tion is opening new insight in the comprehension of this quite complicated 
issue, it is the lack of a theoretical treatment of Compton scattering in the 
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BeppoSAX SVP Observation of Vela. X— 1 
model: constant wabs gausabs gausabs (hlghecut (bknpower) + gaussian) 
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Fig. 3. BeppoSAX Vela X-l br oad-band spectrum fit with three different spectral 
models ( Orlandini et al. , 1998ah . The residual panels show the CRF profile obtained 
by putting the line normalization to zero. See text for the model descriptions. 
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presence of strong magnetic fields the likely culprit of our poor understanding 
of spectral formation XBPs. 



References 

Basko, M., Sunyaev, R., 1976. MNRAS 175, 395. 
Becker, P., Wolff, M., 2005a. ApJ 621, L45. 
Becker, P., Wolff, M., 2005b. ApJ 630, 465. 

Coburn, W., Heindl, W., Rothschild, R., et al., 2002. ApJ 580, 394. 
Filippova, E., Tsygankov, S., Lutovinov, A., et al., 2005. Astron. Lett. 31, 729 
Giacconi, R., Gursky, H., Kellogg, E., et al., 1971. ApJ 167, L67. 
Harding, A., 2003. In: Cusumano, G., Massaro, E., Mineo, T. (Eds.), Pulsars, 

AXPs and SGRs observed with BeppoSAX and other Observatories. Aracne 

Editrice, p. 127. 
Hickox, R., Narayan, R., Kallman, T., 2004. ApJ 614, 896. 
Kreykenbohm, I., Coburn, W., Wilms, J., et al., 2002. A&A 395, 129. 
La Barbera, A., Santangelo, A., Orlandini, M., et al., 2003. A&A 400, 993. 
Makishima, K., Mihara, T., 1992. In: Tanaka, Y., Koyama, K. (Eds.), Frontiers 

of X-ray Astronomy. Universal Academy Press, Tokyo, p. 23. 
Meszaros, P., 1992. High-Energy Radiation from Magnetized Neutron Stars. 

Chicago University Press. 
Meszaros, P., Nagel, W., 1985a. ApJ 298, 147. 
Meszaros, P., Nagel, W., 1985b. ApJ 299, 138. 

Mihara, T., 1995. PhD Thesis, Institute of Physical and Chemical Research, 

Cosmic Radiation Laboratory, Tokyo 
Orlandini, M., 2004. In: Tavani, M., Pellizzoni, A., Vercellone, S. (Eds.), X-ray 

and Gamma-ray Astrophysics of Galactic Sources, Fourth AGILE Science 

Workshop. Aracne Editrice, p. 119 
Orlandini, M., Dal Fiume, D., 2001. In: White, N., Malaguti, G., Palumbo, 

G. (Eds.), X-ray Astronomy '99. AIP Conference Series Vol. 599, p. 283. 
Orlandini, M., Dal Fiume, D., Del Sordo, S., et al., 1999. A&A 349, L9. 
Orlandini, M., Dal Fiume, D., Frontera, F., et al., 1998a. A&A 332, 121. 
Orlandini, M., Dal Fiume, D., Frontera, F., et al., 1998b. ApJ 500, L163. 
Rybicki, G., Lightman, A., 1975. Radiative Processes in Astrophysics. John 

Wiley & Sons. 

Schreirer, E., Levinson, R., Gursky, H., et al., 1972. ApJ 172, L79. 
Shklovskii, I., 1967. ApJ 148, LI. 

Soong, Y., Gruber, D., Peterson, L., et al., 1990. ApJ 348, 641. 
Sunyaev, R., Titarchuk, L., 1980. A&A 86, 121. 

Tanaka, Y., 1986. In: Mihalas, D., Winkler, K. (Eds.), Radiation Hydrody- 
namics in Stars and Compact Objects. Springer, Berlin, p. 198. 
White, N., Swank, J., Holt, S., 1983. ApJ 270, 711. 



8 



